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In this reply to tlie comment by C. R. Willis, we show, by quoting his own statements, that the 
simulations reported in his original work with Boesch [Phys. Rev. B 42, 2290 (1990)] were done for 
kinks with nonzero initial velocity, in contrast to what Willis claims in his comment. We further 
show that his alleged proof, which assumes among other approximations that kinks are initially 
at rest, is not rigorous but an approximation. Moreover, there are other serious misconceptions 
which we discuss in our reply. As a consequence, our result that quasimodes do not exist in the sG 
equation [Phys. Rev. E 62, R60 (2000)] remains true. 
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I. INTRODUCTION 



In his comment fl'], Willis claims that in he used 
a "rigorous" projection operator collective variable for- 
malism for nonlinear Klein-Gordon equations to "prove" 
the continuum sine-Gordon (sG) equation has a long- 
lived quasimode whose frequency ojs = 1.004 is in the 
continuum just above the lower phonon band edge with 
a hfetime (1/ts) = O.OOSFq. In Q we performed two 
numerical investigations which show that neither intrin- 
sic internal modes nor quasimodes exist in contrast to 
previous results. Willis further claims that he "proves" 
that our first numerical investigation could not possibly 
observe the quasimode in principle and our second nu- 
merical investigation actually demonstrates the existence 
of the sG quasimode. He states that his analytic calcu- 
lations and verifying simulations were performed for a 
stationary sine-Gordon soliton at the origin, but that we 
explained our analytic calculations and confirming simu- 
lations in terms of the Doppler shift of the phonons when 
stationary sine-Gordon solitons have a zero Doppler shift. 
In the following we argue that his comment is irrelevant 
because the results Willis refers to in jj] were obtained 
for moving kinks (which we prove by quoting statements 
from '5]), and because his alleged "proof" is not a math- 
ematical proof but an approximate calculation as any 
collective coordinate theory 
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II. ON THE NON-STATIONARITY OF KINKS 

After carefully reading their paper Q several times, 
we are afraid that Willis's comment suffers from severe 
misconceptions and, what is worse, that the author is not 
familiar enough not only with our work j3] but also with 
his own paper 0. Quoting literally from 0, p. 2296, 
first paragraph beginning after Eq. (3.2), it reads: "We 
impose initial conditions in our simulations in two differ- 
ent ways. The first way is by specifying the field [. . . ] 
thus giving the equilibrium kink a nonzero initial veloc- 
ity". And subsequently, three out of four figures of the 
paper refer to this type of simulation. For instance, the 
caption of Fig. 1 reads: "The initial condition for t = 
is a kink with velocity determined by specifying the kink 
shape at. . . ". Therefore, in spite of Willis's claims that 
their simulations were for stationary kinks, that is not the 
case, and he cannot dismiss our explanation of his obser- 
vations based on this. Kinks moved in his simulations 
and correspondingly their phonon spectra were shifted, 
as we explained in Q . 



III. ON WILLIS'S MISCONCEPTIONS 

Given the fact that the author is not aware of his own 
work and results, it is no surprise that he is completely 
unfamiliar with p]. His comment f}\ is full with wrong 
statements about what we reported in that paper. Thus, 
he claims: 

• "What they measured was the phonon absorption 
spectrum of the hnearized sine-Gordon (sG) equa- 
tion. . . " This is simply not correct. We simulated 
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the full sG equation: in page R61 of Q, in the be- 
ginning of section III, this is perfectly clear: "We 
have computed the numerical solution of the per- 
turbed sG equation...". We never linearized any- 
thing and worked with the full equation. 

• "Consequently, their two equations of motion for 
P(t) and l(t) which they state is the 'basis of their 
theoretical analysis' of Ref.[3] have absolutely no 
relevance to our analytic derivation and confirming 
simulations." This is not correct. Actually, those 
equations are the way we find that, were there any 
internal modes or quasimodes in the kink, they 
would be excited parametrically by the ac driv- 
ing. In spite of his claim that he "proved in Sec. Ill 
that an ac driver can not create a sG quasimode" , 
we have shown, by means of those equations, that 
parametric excitation at half the frequency of the 
mode must be observed. We refer the interested 
reader to S learn more about this effect. 
On this issue, it is particularly important to realize 
the inconsistencies in Willis's reasoning: Just af- 
ter Eq. (16) of 1^ Willis claims: "Thus an ac force 
field cannot possibly excite a phonon mode. Con- 
sequently, their first numerical investigation in ref 
[3] could not possibly detect the presence of the 
sG quasimode and thus it has not relevance what- 
soever to the existence or nonexistence of the sG 
quasimode". If this were true (which is not), WiUis 
would be forced to conclude that a biharmonic force 
f{t) = — [ei cos{Ljt) + 62 cos{2ujt + 6)] is also un- 
able to excite the phonons. But this is not the 
case: On the contrary, Willis has recently studied 
0, by using the same approximate method that 
he developed in 1990, the action of this biharmonic 
force and dissipation on the ratchet motion of the 
kink. In this paper, he claims: "We use a rigor- 
ous collective variable for nonlinear Klein-Gordon 
equations to prove that the rectification of the cur- 
rent is due to the excitation of an internal mode 
T{t), which describes the oscillation of the slope of 
the kink, and due to a dressing of the bare kink by 
the ac driver." In Willis's own words, dressing of 
the kinks means phonons excited in the system, in 
blatant contradiction with his remark on our work. 
We also note that Willis's discussion of our theo- 
retical treatment is plagued with his misconception 
about "exact analytical calculations" in his earlier 
paper, which unfortunately are not exact (cf. Sec. 
V below). 

• "The simulations and analysis by the authors in 
Ref. [3] were done for an appreciably discrete sG 
equation". This is not correct. Our discretization 
was of the same order as that of 0: 0.05 in our 
case, 0.02 in theirs [cf. p. 2296, below Eq. (3.1)]. 
If ours is irrelevant, so is the one in 0, and vice- 
versa. The author either did not remember the 
discretization they used, or has not read our paper 



carefully enough to notice ours. 

• On the other hand, Willis does not trust our nu- 
merical simulations and explicitly refers to the "in- 
competent design of their [ours in 3]] simulation", 
claims that all we observe in Fig. 2 of ^ is "a com- 
plicated interference pattern" , arising from "reflec- 
tion of phonons at the boundary at time 200" . We 
were fully aware of the possibility of phonon reflec- 
tion at the boundaries, this being the reason why 
we used free boundary conditions, to minimize this 
effect. Furthermore, we did not trust a visual in- 
spection of our Fig. 2 as he does, and carried out 
a Fourier analysis of the time evolution of the kink 
width, finding only phonons (as Willis points out, 
present in the system due to emission and reflec- 
tion). We did not find any Fourier peak arising 
from a hypothetical "quasimode" , in spite of the 
fact that the author of the comment concedes that 
"the first 200 seconds [actually, time units] of their 
simulation of the width l(t) [actually, the simula- 
tion was for the full sine-Gordon equation, not of 
the width] gives a very good representation of the 
sG quasimode" . How can then be that it does not 
show up in the Fourier spectra, which we report in 
the paragraph below Fig. 2? 



IV. ON FURTHER EVIDENCE ON THE 
NON-EXISTENCE OF QUASIMODES 

Having discussed the lack of familiarity of the author 
of the comment with both his own work and the one he 
comments upon, we would like to refer the reader to re- 
search Willis is not aware of, and that has significantly 
advanced the subject. Indeed, one of us (NRQ) has pub- 
lished a very relevant investigation on this subject, 0, 
showing that the spectrum mentioned in the comment, 
see Eq. (15), and which we used in our own paper [Eq. 
(15) of Q]] is not correct, and in turn the parenthesis 
should read {n — 1)tt/L for small n/L. This paper is 
very pertinent to this discussion because the author of 
the comment seems to be very worried that we are not 
aware of the "phonon dressing" of the kink. This work of 
Quintero and Kevrekidis is the best way to understand 
what are the consequences of the phonons on the kink 
behavior. 

Furthermore, there was only one other paper in the lit- 
erature that ever observed something similar to a "quasi- 
mode" |lO|[ . in fact using simulations of kinks driven by 
constant forces. This only independent confirmation of 
the "quasimode" was also proven wrong by us in (lli ]. 
In this case, it is not a matter of simulations and how 
they measured the "quasimode" : In the above reference 
we showed that the analytical calculations in this other 
paper were not correct. We encourage the reader to take 
this further report on the absence of quasimodes into ac- 
count. 
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V. ON THE APPROXIMATE CHARACTER OF 
COLLECTIVE COORDINATE THEORIES 

Finally, Willis does not seem to be informed about the 
nature and rigor of collective coordinate calculations. He 
states in the abstract that "... prove the continuum sine- 
Gordon equation has a long lived quasimode..." . Unfor- 
tunately and contrary to his belief, collective coordinate 
calculations can not prove anything in the mathemati- 
cal sense, in so far as they are approximate. Even if 
he tries to keep his calculation (which, by the way, he 
reproduces literally from 0] , therefore adding no new ev- 
idence at all in favor of his claims) exact, he cannot, and 
he has to make a number of approximations including 
two linearizations and assuming X and its derivative to 
be exactly zero (see p. 2293 of 0, 2nd paragraph on 
the right). Willis goes as far as claiming: "their Eq. (3) 
for l{t) is incorrect because it contains none of the many 
terms proportional to x(t) that appear in the exact equa- 
tion of motion. . . " Our equation is completely correct, 
consistent with our derivation, and confirmed by numer- 
ical simulations. Our result is different of his simply be- 
cause neither of them is exact. They begin with different 
Ansdtze and subsequently lead to different approximate 
equations. Collective coordinate and related calculations 
are interesting, often times very accurate and useful, but 
by no means exact and, mathematically speaking, con- 
stitute no proof. Indeed, it has been shown that such a 
procedure may lead to wrong predictions: For a detailed 
discussion, see jj]. 

In this respect, we want to stress that when Willis 
uses the paper |l2| in favor of his argument, it is again 
an approximate result and a conjecture. It is true that 
Kalbermann calculates some solutions analytically, but 
he has to resort to numerical simulations to make his 
point. In addition, the relationship of Kalbermann's 
"wobble solution" with Willis's "quasimode" is dubious 
at best: Kalbermann finds a continuum of such possible 
quasimodes, which is not surprising because his wobble 
solution is a combination of a kink and a breather, and 



therefore is a one-parameter family of solutions; a longer 
discussion of this work is out of the scope of this paper, 
but we want to stress that he never finds a specific quasi- 
mode as Willis claims, and Kalbermann's only example 
has a frequency below the phonon band, very different 
from Willis's quasimode. Furthermore, this wobble is, in 
Kalbermann's own words, only "apparently stable" and 
all he can be positive about is that "the kink appears to 
be decaying to a wobble" . Therefore, there is no such 
thing as a mathematical proof of the existence of "quasi- 
modes" in this paper either. 



VI. CONCLUSION 

Based on all the arguments summarized above, we are 
afraid that the preceding comment T] is a combination 
of wrong memories and misconceptions about our work. 
Willis is not aware that his own work did not deal with 
stationary kinks, he is not familiar with the fact that the 
simulations in |2| are as discrete as ours, does not un- 
derstand that collective coordinate theories are not rig- 
orous mathematical proofs, is not aware of further devel- 
opments on this question, and reads our paper Q in an 
unscientifically biased manner. Therefore, Willis's com- 
ment does not change our conclusions, and our result that 
quasimodes do not exist in the sG equation Q stands in 
its own right. 
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